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ABSTRACT 

The  rate  at  which  the  solid  phase  develops  in  a  fluid,  moving 
with  Blasius  flow  over  a  chilled  flat  plate  is  analyzed.   The  time 
dependent,  local,  solid  layer  thickness  is  determined  as  a  function 
of  the  pertinent  physical  properties s  the  plate  surface  temperature 
and  a  one  spatial  dimension  heat  flux.   Experimental  equipment  to 
measure  the  heat  flux,  transient  ice  growth  and  steady  state  ice 
profile  is  described*   The  spatial  variation  of  the  steady  state  ice 
profile  is  observed  with  respect  to  the  streamwise  coordinate. 
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NOMENCLATURE 


a1  dimensionless  Blasius  stream  function 

C.  constant  of  integration 

h  latent  heat  of  fusion 

Li 

k  thermal  conductivity  of  the  liquid  phase 

Li 

k  thermal  conductivity  of  the  solid  phase 

P  Prandtl  Number 

q  heat  flux  of  the  liquid  phase  in  the  y  direction 

q  heat  flux  of  the  solid  phase  in  the  y  direction 

R  Reynolds  Number 
e 

t  time 

t  parameter  with  the  dimension  of  time 
s 

Tf  fusion  temperature  of  the  solid  phase 

T  local  temperature  of  the  liquid  phase 

Li 

T  local  temperature  of  the  solid  phase 

T  temperature  of  the  cold  plate  surface 
w 

Tgo  ambient  liquid  phase  temperature 

Uoo  free  stream  velocity  in  the  x  direction 

x  dimensional  position  coordinate  tangent  to  the 
plate  in  the  direction  of  flow 

y  dimensional  position  coordinate  normal  to  the  cold 
plate  horizontal  surface 


Greek  Symbols 

TO  transient  solid  phase  thickness 

Vs  steady  state  solid  phase  thickness 

7>  kinematic  viscosity 

J  dimensionless  thickness 

J*L  mass  density  of  the  liquid  phase 

t  dimensionless  time 
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1.    Introduction 

In  recent  years  there  has  been  increased  attention  directed  toward 
solving  the  problems  of  heat  transfer  through  a  moving  fluid  undergoing 
a  phase  change.   The  greater  part  of  the  work  in  this  area  is  concen- 
trated in  the  area  of  boiling;  its  various  modes,  bubble  behavior  and 
interactions.   On  the  other  hand,  the  rate  of  solid  phase  development 
in  a  fluid  and  the  geometry  of  the  steady  state  profile  have  been  sub- 
jected to  some  new  investigations.   Much  of  the  mounting  interest  in 
"phase -change"  problems  can  be  found  in  the  aerospace  and  nuclear  power 
industries.   In  the  analysis  of  airfoils,  consideration  must  be  given 
to  ice  build-up.   Proper  design  of  space  vehicles  demands  hydraulic 
systems  and  heat  exchangers  that  are  free  of  unwanted  phase  changes. 
Nuclear  reactor  engineers  are  studying  the  complexities  associated 
with  using  liquid  metals  in  coolant  loops.   Heat  transfer  with  a  phase 
change  is,  therefore,  one  of  the  problems  that  has  become  more  acute 
as  modern  technology  advances. 

One  of  the  earliest  treatments  of  the  problem  of  heat  conduction 
through  a  solid  in  the  presence  of  a  change  of  phase  was  Stefan's  [1]* 
classical  paper  on  the  formation  of  polar  ice.   He  determined  the  posi- 
tion of  the  freezing  line  by  formulating  a  one  spatial  dimension,  heat 
conduction  problem  for  a  semi-infinite  solid.   Stefan  assumed  the  water 
mass  was  at  its  solidification  temperature  and  he  obtained  an  equation 
for  the  position  of  the  interface.   This  equation  was,  however,  in  a 


* 
Numbers  in  brackets  refer  to  references  listed  in  the  Bibliography 
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complex  mathematical  form,   A  simplified  version  of  that  equation  developed 
when  the  heat  capacity  of  the  ice  was  neglected,   Stefan's  approximate 
solution  to  the  original  equation  stipulated  a  linear  decrease  in  the 
solid  phase  temperature  between  the  melting-point  temperature  at  the 
interface  and  an  assumed  solid  phase  temperature  of  zero  at  the  water- 
line  of  the  ice  mass.   Neumann  [2]  provides  an  exact  solution  to  the 
Stefan  problem  without  restricting  the  liquid  mass  temperature  to  a 
constant  or  having  to  assume  a  temperature  gradient  in  the  solid  phase. 
Even  in  a  non-flow  situation,  the  phase-change  heat  transfer  problem 
becomes  non-linear  and  requires  special  solutions  that  cannot  be  super- 
posed.  The  "melting-freezing'"  problem  in  a  fluid  flow  field  is  diffi- 
cult to  analyze  because  the  spatial  position  and  profile  of  the  phase 
change  boundary  are  not  known  a  priori  and  therefore,  must  be  found  as 
part  of  the  problem  solution, 

A  review  of  recent  literature  discloses  that  mathematical  analysis 
using  approximate  methods  is  employed  in  the  works  of  Goodman  [3,4,5], 
Biot  [6,7],  Libby  and  Chen  [8],  and  Murray  and  Landis  [9].   Goodman 
and  Shea  [3]  presented  an  analysis  of  the  melting  problem  in  a  finite 
slab.   Cho  and  Sunderland  l"10]  used  Goodman's  integral  technique  to 
find  the  temperature  distribution  in  their  work  on  melting  finite  slabs. 
Tien  and  Yen  [11]  applied  Goodman's  methods  to  obtain  an  approximate 
solution  to  a  Stefan  problem  in  which  a  semi-infinite  slab  melts  at  its 
lower  face  and  the  liquid  phase  is  subjected  to  a  negative  temperature 
gradient.   Assuming  the  liquid  density  is  inversely  proportional  to  the 
temperature,  the  system  is  seen  to  be  inherently  unstable  and  the  heat 
transfer  in  the  liquid  is  enhanced  by  natural  convection. 
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Liquid  solidification  in  a  circular  tube  under  steady  state  condi- 
tions was  analyzed  by  Zerkle  and  Sunderland  [12]  and  equations  for  the 
temperature,  velocity  and  pressure  distributions  in  a  liquid  flowing 
through  a  circular  tube  with  solidification  at  the  inner  wall  surface 
were  derived.   In  the  analysis  it  was  assumed  free  convection  within 
the  liquid  was  negligible.   The  experimental  data  presented  showed  a 
significant  variation  between  theoretical  predictions  and  experimental 
results  when  convection  exists.   Transient  freezing  of  a  liquid  flowing 
inside  a  circular  tube  was  investigated  analytically  by  Ozisik  and 
Mulligan  T13].   They  assumed  the  liquid  flow  characteristics  inside  the 
tube  could  be  represented  by  a  slug-flow  velocity  profile.   Their  other 

assumptions  were  similar  to  those  of  Zerkle  and  Sunderland  T12], 

ii 

Ozisik  and  Mulligan  obtained  closer  correlation  to  Zerkle  and  Sunder- 
land's experimental  data  for  steady-state  solid-liquid  profiles  and 
total  steady-state  heat  transfer  rates  than  did  Zerkle  and  Sunderland. 
The  results  of  the  transient  part  of  Ozisik  and  Mulligan  s  analysis  was 
not  correlated  with  experimental  data  because  no  published  data  existed 
for  their  situation.   They  limited  the  applicability  of  their  results  to 
small  rates  of  change  of  thickness  of  the  solid  layer  with  time  and  with 
distance  along  the  tube. 

The  freezing  of  fluids  in  forced  flow  over  a  flat  plate  has  been 
analyzed  mathematically  by  Beaubouef  and  Chapman  [14].   Their  exact 
numerical  solution  for  solid  phase  thickness  as  a  function  of  time  in- 
cludes the  heat  capacity  of  the  ice  layer.   It  compares  closely  with 
the  approximate  solution  of  Lapadula  and  Mueller  F 15 ]  who  used  Biot's 
methods . 
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The  unsteady  formation  of  ice  in  a  plane  stagnation  flow  was  studied 
by  Yang  [16].   The  problem  dealt  with  an  instantaneously  cooled  plate 
normal  to  an  isothermal  flow  of  liquid.   Short  and  long  time  solutions 
were  obtained  analytically  for  the  ice  thickness  profile.   Velocity 
profile  distortions  caused  by  the  advance  of  the  ice  front  were  also 
discussed. 

Published  experimental  results  relevant  to  solidification  on  a 
flat  surface  geometry  are  somewhat  limited.   The  work  of  Savino  and 
Siegel  r 17 ,  18]  provides  much  of  the  available  experimental  information. 
They  did  not  consider  the  details  of  the  fluid  flow  field  in  their 
analysis  and  correlations  with  experimental  results.   The  heat  capa- 
cities of  the  plate  and  the  solidified  material  were  accounted  for  in 
reference  18. 

This  thesis,  therefore,  was  undertaken  to  compare  the  experimental 
findings  of  Siegel  and  Savino  using  a'plate  of  similar  geometry,  but 
larger  in  size.   Instrumentation  was  provided  to  gain  information  on 
transient  solidification  in  the  flow  field  as  a  function  of  time  and 
also  as  a  function  of  the  streamwise  spatial  coordinate.   The  capability 
of  direct  measurement  of  the  thickness  contour  of  the  steady  state  ice 
profile  was  included  in  the  experimental  equipment  design.   Further,  a 
different  analysis  was  developed  for  predicting  the  transient  thick- 
ness of  the  solidified  layer. 
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2.    Theoretical  Analysis 

The  problem  being  considered  is  the  freezing  of  a  fluid  in  steady, 
laminar,  forced  flow  over  a  flat  plate  that  has  been  suddenly  cooled. 

The  basic  assumptions  used  are  similar  to  those  stated  in  refer- 
ences 14,  17  and  18.   They  are: 

a.  The  thickness  of  the  deposited  layer,  "W  ,  is  such  that 
heat  conduction  in  the  solid  layer  is  considered  one-dimensional. 

b.  Blasius  flow  is  assumed  so  that  the  fluid  boundary  layer 
is  not  affected  by  the  solid  phase  geometry. 

c.  The  temperature  of  the  cold  plate,  T  ,  and  the  solid- 
fluid  interface  temperature,  Tf,  are  considered  constant. 

d.  The  physical  properties  of  both  the  liquid  and  solid 
phases  are  constant. 

e.  The  heat  capacity  of  the  ice  layer  is  negligible. 

The  one-dimensional  model  and  some  of  the  notation  employed  are 
shown  in  Figure  1.   A  warm  liquid,  moving  with  Blasius  flow,  passes 
over  a  flat  plate.   The  plate  is  chilled  and  freezing  is  assumed  to  be- 
gin (t=0)  as  the  plate's  upper  surface  reaches  the  temperature  of  solid- 
ification of  the  liquid.   During  the  freezing  process  the  latent  heat 
of  fusion,  h  ,  is  released  at  the  liquid-solid  interface.   The  solid 

Li 

phase  develops  until  it  reaches  a  steady-state  thickness,  Vc    •      It  is 
assumed  that  this  thickness  will  never  become  so  great  as  to  disturb 
the  Blasius  flow. 

Neglecting  the  heat  capacity  of  the  ice,  the  one  -  dimensional  heat 
flux  for  the  solid  phase  is: 
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where  the  solid  phase  thickness  f)     ,  is  a  function  of  time  and  the 
streamwise  spatial  coordinate,  x. 

An  interface  energy  balance  yields: 

For  Blasius  flow  the  temperature  gradient  is  given  by  reference 
19  as: 


where  Q1   =  0.33  Z  PR'/3 

Substituting  equation  (3)  into  equation  (2)  and  solving  for  the 
freezing  rate 

For  steady  state  conditions 

Then  equation  (4)  can  be  solved  for  the  steady  state  thickness 
rj  -Ks      Of  ""Tvvi^ 


K        (to-TJ  ^?r.332)Pt* 

Substituting  in  the  Reynold's  number 


get 


(5) 
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Equation  (4)  can  then  be  transformed  into  non-dimensional  form  by 


dividing  by  Vs     : 


K    Re*  (532)PR* 


y\itj    at\y   ?L hL  \  %*■/  ?      ?LhL  x    $ 


(t~-t<) 


Rs  I  if  ~  iw 


AbeV   $ 


2. 


&  _  &.  (T~-Tf)(.55Z)fl£PeS  n  1 


Tf-Tw 


fc 


-'] 


(6) 


Defining  a  parameter  whose  dimension  is  time  as 


.  _  Ak/JL_ 


non-dimensional  time  becomes 


and  non-dimensional  thickness  is 

Therefore,  equation  (6)  becomes 

*  (i)-4-i 


at 


l-I 


or 


=  dr 


Integrating  equation  (7),  we  get 


(7) 


(8) 
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To  evaluate  the  constant  of  integration  the  following  initial 
condition  is  used: 


at 


r=o,  M  =  o 


Then  from  equation  (8)  get 

C,-l 

and 


or 


Figure  2  is  a  graph  of  equation  (9) . 

The  above  solution  can  be  shown  to  be  the  limiting  case  for  a  solid 
phase  heat  capacity  equal  to  zero  in  the  analyses  of  Beaubouef  and 
Chapman  T14]  and  Savino  and  Siegel  [17]. 

This  analysis  is  unique  in  its  straightforward  derivation  of  the 
closed  form  solution  as  compared  with  the  equations  given  in  references 
14  and  17.   The  above  solution  also  has  a  distinct  advantage  in  its  ease 
of  computation. 
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3.    Experimental  Equipment  and  Procedure 

The  experimental  equipment  shown  schematically  in  Figure  3 
utilized  a  variable  speed,  axial  flow  pump  to  control  the  flow  of 
ordinary  tap  water  over  the  test  plate.   The  pump  took  suction  from 
a  reservoir.   Re-circulation  of  a  large  water  volume  during  each  test 
run  did  not  change  the  temperature  of  the  liquid  because  of  the  large 
thermal  capacity  of  the  water  mass.   From  the  pump,  the  water  passed 
through  an  A.S.M.E.  calibrated  orifice,  a  diffuser  and  a  series  of  flow 
straighteners  consisting  of  a  metal  honeycomb  and  wire  screens.   The 
test  section  of  the  water  channel  was  nominally  one  foot  square  in 
cross-section.   From  the  channel,  the  liquid  flowed  into  a  back-water 
tank  which  was  used  as  an  aid  to  flow  control.   The  water  passed  from 
the  tank,  through  a  control  valve  and  a  rotometer,  back  to  the  reservoir. 
The  flow  straighteners  in  the  channel  upstream  of  the  test  section  pro- 
vided a  laminar,  uniform  water  velocity  profile  across  the  plate  sur- 
face.  The  test  section  itself  was  a  clear  plastic  rectangular  chan- 
nel (see  Figure  4).   The  support  frame  for  the  test  section  was  adjust- 
able so  as  to  provide  a  zero  pressure  gradient  in  the  flow  field  along 
the  bottom  of  the  channel. 

The  refrigerant  was  re-circulated  through  a  closed  loop  by  a  one- 
half  horsepower  commercial  compressor  unit.   The  boiling  chamber  was 
soldered  to  the  bottom  of  the  test  plate.   Freon  flow  was  regulated 
through  a  manually  controlled  expansion  valve.   The  temperature  of  the 
refrigerant  was  adjusted  in  each  test  run  by  control  of  the  pressure 
in  the  boiling  chamber. 

The  plate  was  machined  from  oxygen  free,  high  capacity,  pure  copper 
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12  inches  long,  6  inches  wide  and  2  inches  thick.   Copper  was  selected 
as  the  test  material  because  of  its  extremely  high  thermal  conductivity. 
This  property,  acting  in  conjunction  with  the  plate  thickness  was  deemed 
sufficient  to  insure  a  constant  plate  temperature  at  its  top  surface. 
The  upper  surface  of  the  test  plate  was  machined  flat  and  highly  polish- 
ed to  minimize  the  possibility  of  any  flow  disturbances  across  it. 

The  test  plate  and  refrigerant  boiling  chamber  were  supported  by  a 
frame  designed  to  reduce  conduction  heat  losses  to  a  minimum  and  to  pro- 
vide position  adjustment  such  that  the  top  surface  of  the  plate  lay  in 
the  plane  of  the  channel  bottom  (see  Figure  4).  A  rubber-base,  water- 
proof sealant  of  low  thermal  conductivity  was  used  to  fill  the  gap  be- 
tween the  plate  edge  and  the  plastic  of  the  channel  bottom.  All  other 
exposed  sides  of  the  plate  and  the  boiling  chamber  were  covered  with  a 
styrofoam  insulating  material. 

Instrumentation  of  the  test  plate  consisted  of  nine  stainless  steel 
sheathed,  ceramic  insulated,  ungrounded,  copper-constantan  thermocouples, 
The  outside  diameter  of  the  sheathing  was  1/16  inch  and  the  thermocouple 
wires  were  number  30  AWG,  with  a  nominal  diameter  of  10  mills.   The 
location  of  the  thermocouples  is  indicated  in  Figure  5.   This  particular 
arrangement  provided  a  temperature  profile  across  the  central  area  of 
the  top  horizontal  surface  of  the  plate.   Five  evenly  spaced  points 
across  the  thickness  of  the  plate  near  the  vertical  centerline  of  the 
mass  of  the  plate  were  selected  to  provide  heat  flux  information.   These 
thermocouple  positions  were  purposely  staggered  slightly  off  the  plate 
centerlines  to  reduce  mutual  interference  and  minimize  discontinuities 
in  the  heat  flow. 
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Benefiting  from  the  experience  of  Savino  and  Siegel  [17],  a  thermo- 
couple probe  (Figure  6)  was  used  to  measure  the  ice  layer  thickness. 
Its  construction  was  similar  to  that  described  in  reference  17.   The 
probe  consisted  of  butt  welded  copper-constantan  wires  of  .003  inch 
diameter.   The  wires  were  supported  by  two  sewing  needles  mounted  at  90° 
angles  to  the  bottom  of  a  £-inch  stainless  steel  rod.   The  needle  moun- 
ting was  designed  to  reduce  flow  disturbance  around  the  thermocouple 
junction.   The  rod  extended  up  to  the  top  of  the  channel  where  it  was 
attached  to  a  micrometer  head  through  a  sleeve  and  cross-bar  device. 
The  micrometer  spindle  was  mounted  in  a  ball  bearing  attached  to  the 
cross  channel  carriage,  through  which  the  sleeve  and  probe  traveled. 
The  entire  device  enabled  the  micrometer  head  to  raise  and  lower  the 
thermocouple  while  maintaining  a  preset  alignment  of  the  needles  with 
respect  to  the  fluid  flow.   The  micrometer  head  was  attached  to  a  two 
dimensional  traversing  carriage.   The  carriage  frame  was  mounted  on  the 
test  section  with  leveling  screws.   Thus,  the  thermocouple  junction 
enjoyed  movement  in  three  degrees  of  freedom  and  could  be  accurately 
positioned  at  any  specified  height  above  the  plate  as  well  as  at  any 
position  over  the  horizontal  area  of  the  test  plate.   Machinist  scales 
and  indicators  were  attached  to  the  traversing  carriage  to  position  the 
thermocouple  at  various  pre-selected  points  on  the  plate's  horizontal 
surface  during  successive  test  runs.   As  Savino  and  Siegel  T17]  noted: 
The  thermocouple  and  its  needle  supports  were  sufficiently  small  so  that 
no  local  distortions  of  the  interface  were  observed  when  the  thermo- 
couple was  placed  in  contact  with  it. 

The  reference  junction  for  all  thermocouples  was  boiling  nitrogen 
at  atmospheric  pressure.   All  test  data  was  observed  away  from  the  edges 
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of  the  plate  where  the  heat  transfer  was  more  closely  one-dimensional 
(a  basic  assumption  in  the  theoretical  analysis).   A  procedure  described 
by  Savino  and  Siegel  r 17 ]  was  used  to  measure  instantaneous  thickness 
during  a  test  run.   The  thermocouple  junction  was  set  0.005  inch  above 
the  plate,  and  the  probe  temperature  was  monitored  using  a  taoe  output 
recorder.   As  the  ice  layer  developed,  the  probe  temperature  decreased 
and  at  32°F.  the  ice  layer  contacted  the  probe.   The  thermocouple  junc- 
tion was  then  raised  0.005  inch  rapidly  so  that  it  would  not  freeze  in 
the  ice.   The  recorder  then  registered  a  step  increase  in  the  fluid 
temperature  and  the  cycle  was  repeated.   The  recorder  output  indicated 
a  series  of  steps  and  the  chart  recorded  the  time  increment  correspond- 
ing to  each  0.005  inch  of  growth  of  the  ice  layer. 

Direct  measurement  of  the  thickness  of  the  ice  contour  was  under- 
taken through  use  of  the  two  dimensional  traversing  carriage.   A  length 
of  drill  rod  was  substituted  for  the.  thermocouple  probe.   The  lower  end 
of  the  rod  was  shaped  into  a  semi-sphere.   The  machinist  scales  and  in- 
dicators, mounted  on  the  carriage  frame  in  the  direction  of  flow  and 
normal  to  it,  provided  a  means  to  position  the  depth  indicating  rod 
over  any  predetermined  plate  position.   The  micrometer  head  was  set  at 
zero  when  the  rod  rested  on  the  plate  surface.   In  this  configuration 
the  traversing  carriage  indicated  the  ice  thickness  at  selected  uni- 
form intervals  over  the  horizontal  surface  of  the  layer. 
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4.    Discussion  of  Theoretical  and  Experimental  Results 

The  theoretical  analysis  assumed  that,  at  zero  time,  the  plate  sur- 
face over  which  the  water  was  flowing  was  at  the  freezing  temperature 
and  ice  began  to  form.   Under  experimental  conditions,  ice  formed  rapidly 
and  evenly  over  the  horizontal  surface  of  the  plate.   Five  seconds 
elapsed  between  formation  of  the  initial  patch  of  ice  on  the  upstream  edge 
until  the  plate  was  completely  covered.   Ice  formed  along  the  direction 
of  flow.   This  fact,  coupled  with  the  observation  that  the  refrigerant 
entered  the  boiling  chamber  at  the  upstream  edge  of  the  plate,  makes  ap- 
parent the  existence  of  a  layer  of  subcooled  water  in  contact  with  the 
plate  prior  to  the  initial  ice  nucleation.   Positive  verification  of  this 
phenomenon  was  not  obtained  because  the  thermocouple  recorders  available 
for  use,  at  the  time  of  the  experiment,  did  not  possess  the  required 
sensitivity.   Electromotive  forces  of  the  order  of  5  millivolts  were  en- 
countered when  liquid  nitrogen  was  used  as  a  thermocouple  reference  junc- 
tion.  A  full  scale  value  of  5  millivolts  on  a  recorder  strip  chart  re- 
sulted in  poor  temperature  discrimination  near  the  freezing  level.   Trans- 
ient temperature  data  from  the  plate  and  the  probe  did  not  provide  a 
valid  comparison  for  the  theoretical  analysis  because  of  the  same  re- 
corder limitations. 

Meaningful  data  was  hampered  by  particles  of  rust  and  debris  that 
deposited  on  the  bottom  of  the  test  section.   The  sources  of  this  foreign 
matter  and  the  oil  droplets  in  the  water  was  traced  to  the  piping  sys- 
tem and  to  leakage  of  lubricating  oil  from  the  axial  flow  pump  respec- 
tively. 

The  steady  state  ice  layer  was  nominally  .1  inch  thick  over  the 
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surface  of  the  plate.   Edge  effects  were  confined  to  a  1/8  inch  peri- 
meter around  the  downstream  edge  and  the  two  edges  parallel  to  the  flow. 
Streaks  in  the  ice  layer  extended  from  the  upstream  edge  for  \   inch  in 
the  flow  direction.   The  streaks  were  caused  by  small  irregularities  in 
the  sealant  applied  at  the  edge  and  to  debris  on  the  bottom  of  the  test 
section.   Localized  depressions  in  the  ice  layer  thickness  occurred 
around  the  sheathed  thermocouples  imbedded  in  the  plate's  surface. 
These  irregularities  became  less  pronounced  as  the  ice  thickness  in- 
creased.  In  all  other  areas  of  the  plate  the  ice  was  clear  and  its 
surface  was  very  smooth. 

While  measuring  the  thickness  of  the  steady  state  ice  layer,  the 
rod  was  observed  to  cause  depressions  in  the  ice  layer.   If  the  rod  was 
left  in  the  proximity  of  the  layer  it  would  distort  the  flow  to  the  ex- 
tent that  ridges  would  develop  on  the  surface  of  the  ice. 

Much  time  was  expended  in  obtaining  the  optimum  refrigerant 
charge  for  the  system.   Changes  in  flow  velocity  had  great  effect  on 
the  cooling  rate  produced  by  the  freon. 


24 


5.    Conclusions 

The  object  of  the  theoretical  analysis  was  to  develop  a  means  for 
predicting  the  transient  growth  of  the  solid  phase  of  a  liquid  flowing 
over  a  chilled  flat  plate.   The  object  was  achieved  by  making  two  basic 
assumptions:   (1)  the  liquid  moved  over  the  plate  with  Blasius  flow  and 
(2)  conduction  in  the  solid  was  one  dimensional  and  the  heat  capacity  of 
the  ice  was  neglected. 

The  layer  of  subcooled  water  in  contact  with  the  plate  prior  to 
ice  nucleation  must  be  accounted  for  to  make  experimental  data  corres- 
pond with  the  analysis. 

The  recording  equipment,  in  its  present  configuration,  will  not 
provide  usable  data  on  transient  temperatures  in  the  plate.   The  thermo- 
couple probe  record  did  not  provide  an  accurate  indication  of  solid 
phase  development.   The  rod  used  to  measure  the  steady  state  ice  layer 
thickness  distorted  the  flow  and  it  conducted  heat  toward  the  ice  sur- 
face. 

An  acceptable  steady  state  ice  layer  thickness  can  be  obtained  by 
increasing  the  capacity  of  the  refrigerant  boiling  chamber.   In  its 
present  design,  almost  half  of  the  area  of  the  chamber  that  is  in  con- 
tact with  the  plate  is  used  for  baffles  and  flow  straighteners .   More 
flexibility  for  the  refrigerant  system  to  meet  a  wide  spectrum  of 
cooling  rates  would  be  a  desirable  feature  in  experimental  equipment 
of  this  kind. 

The  difference  in  conductivity  of  the  bulky  thermocouples  im- 
bedded in  the  surface  of  the  plate  and  the  solder  that  surrounds  them 
cause  noticeable  distortions  in  the  temperature  distribution  and  the 
ice  layer. 
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6.    Recommendations  for  Future  Work 

Install  a  filter  bed  between  the  pump  and  the  test  section  to  remove 
the  suspended  dirt,  rust  particles  and  oil  seepage  from  the  pump  lubri- 
cating system. 

Obtain  Blasius  flow  over  a  greater  range  of  water  channel  veloci- 
ties by  modifications  in  the  design  and  combinations  of  the  metal  honey- 
comb and  wire  screens  used  to  straighten  the  flow. 

A  table  of  test  section  flow  velocities  versus  pump  speed,  pump 
control  board  settings  and  depth  of  water  in  the  channel  would  aid 
investigators  who  plan  to  use  the  channel  for  further  experiments. 

Use  a  multi-channel  thermocouple  strip  recorder  or  a  digital 
voltmeter  coupled  with  an  automatic  sequential  scanning  switch  to 
obtain  transient  temperatures. 

Investigate  the  feasibility  of  using  other  than  boiling  liquid 
nitrogen  as  a  reference  junction.   For  example,  an  ice  bath  may  work 
with  a  digital  voltmeter  or  liquid  oxygen  may  provide  a  more  stable 
reference. 

Increase  the  volume  of  the  refrigerant  boiling  chamber  and  reduce 
the  thickness  of  the  baffles.   Install  a  receiver  in  the  refrigeration 
unit  to  increase  the  volume  of  freon  in  the  system.   This  in  turn 
would  improve  the  flexibility  of  the  system  to  respond  to  variations 
in  plate  cooling  rates  imposed  by  an  investigator. 

Reduce  the  size  of  the  thermocouples  on  the  surface  of  the  plate. 
The  difference  in  conductivity  of  the  sheathing  and  copper  affect  the 
shape  of  the  ice  layer.   Small  diameter  wire  thermocouples  would  re- 
quire more  shallow  grooves  and  less  solder.   They  have  the  added 
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advantage  of  a  faster  time  response  than  the  sheathed  type. 

Insulate  the  tip  of  the  rod  used  for  steady  state  profile  measure- 
ment, with  a  material  of  low  thermal  conductivity.  Streamline  both  the 
rod  and  the  probe  shafts  to  reduce  the  flow  disturbances. 

Obtain  experimental  data  holding  one  of  the  following  variables 
constant  while  measuring  the  other  two: 

a.  water  flow  rate 

b.  refrigerant  flow  rate 

c.  position  on  the  plate  horizontal  surface. 
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-   FIGURE   A.      TEST   SECTION 
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NOTE:  ALL  THERMOCOUPLE  LOCATIONS  AND  PLATE  DIMENSIONS 
ARE  IN  INCHES, 

FIGURE  5.  TEST  PLATE  INSTRUMENTATION 
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